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ABSTRACT 

Context. Strong lens candidates have been newly identified within the COSMOS field. We present VLT/FORS 1 spectroscopic follow- 
up observations and HST/WFPC2 imaging of the system COSMOS 5921+0638, which exhibits quadruply lensed images and a perfect 
Einstein ring. 

Aims. We investigate the nature of COSMOS 5921+0638 by studying its photometric, spectroscopic and physical properties. 
Methods. By analyzing our VLT/FORS 1 spectroscopy and Subaru/CFHT/HST imaging of COSMOS 5921+0638, we completed both 
an environmental analysis and detailed analytical and grid-based mass modeling to determine it properties. 

Results. We measured the redshifts of the lensing galaxy in COSMOS 5921+0638 (z ; =0.551±0.001) and 9 additional galaxies in 
the field (5 of them at z~0.35). The redshift of the lensed source was inferred by identifying a candidate Lyo- line at z s =3. 14+0.05. 
The mass modeling reveals the requirement of a small external shear (y=0.038), which is slightly larger than the lensing contribution 
expected by galaxy groups along the line-of-sight obtained from the zCOSMOS optical group catalog (/c ?roapj ~0.01 and y groups ~0.00S). 
The estimated time-delays between the different images are of the order of hours to half a week and the total magnification of the 
background source is yui»150. The measured mass-to-light ratio of the lensing galaxy within the Einstein ring is M/L B «:8.5±1.6. 
Anomalies are observed between the measured and expected flux ratios of the images of the background AGN. 
Conclusions. Our analysis indicates that the ring and point-like structures in COSMOS 5921+0638 consist of a lensed high redshift 
galaxy hosting a low luminosity AGN (LLAGN). The observed flux ratio anomalies are probably due to microlensing by stars in 
the lensing galaxy and/or a combination of static phenomena. Multi-epoch, multi-band space-based observations would allow us to 
differentiate between the possible causes of these anomalies, since static and/or dynamic variations could be identified. Because of its 
short time-delays and the possibility of microlensing, COSMOS 5921+0638 is a promising laboratory for future studies of LLAGNs. 

Key words, galaxies: individual: COSMOS 5921+0638 - galaxies: quasars: individual: COSMOS 5921+0638 - gravitational lensing 
- cosmology: observations 



1. Introduction 



* Based on observations made with ESO Telescopes at Paranal 
Observatory programme IDs: 077A-0473(A) and 175.A-0839(B,D) 
and with NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Institute. STScI is operated by the as- 
sociation of Universities for Research in Astronomy, Inc. under the 
NASA contract NAS 5-26555; also based on data collected at: the 
Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan; the European Southern Observatory under 
Large Program 175.A-0839, Chile; Kitt Peak National Observatory, 
Cerro Tololo Inter-American Observatory, and the National Optical 
Astronomy Observatory, which are operated by the Association of 
Universities for Research in Astronomy, Inc. (AURA) under coopera- 
tive agreement with the National Science Foundation; and the Canada- 
France-Hawaii Telescope operated by the National Research Council of 
Canada, the Centre National de la Recherche Scientifique de France and 
the University of Hawaii. 



Gravitationally lensed quasars are powerful observational tools 
in cosmology. The intrinsic variability of AGNs, provides a way 
to measure the time-delays between multiple quasar images (an 
updated list of lensed quasars' time-delays was established by 
Oguri 2007). By coupling these measurements with theoreti- 
cal models for the lensing potential, the Hubble constant (Ho) 
can be derived (Refsdal 1964). They allow us to investigate 
the amount and distribution of matter in lensing galaxies re- 
gardless of whether the matter is luminous or not. This prop- 
erty ensures that lensed quasars are ideal laboratories for prob- 
ing for the dark matter distribution and dark matter substruc- 
tures in and around lensing galaxies (Mao & Schneider 1998; 
Chiba 2002; Yonehara et al. 2003; Kochanek & Dalai 2004; 
Keeton & Moustakas 2008; Trott et al. 2008). Furthermore, the 
natural magnification of background sources offers a means of 
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studying the properties of objects in the high redshift universe. 
In particular, if lensed AGN images are affected by brightness 
fluctuations produced by individual stars in lensing galaxies (mi- 
crolensing), the inner structure (i.e., milliparsec scales) of the 
AGN accretion disks can be probed (Wambsganss et al. 1990; 
Yonehara et al. 1999; Kochanek 2004; Anguita et al. 2008b; 
Eigenbrod et al. 2008). 

Strong lens samples (e.g., Myers & The Class Collaboration 
2001; Fassnacht et al. 2004; Cabanacetal. 2007; Bolton etal. 
2008; Inada et al. 2008; Faure et al. 2008) allow us to study in 
a statistical way the properties of lensing galaxies (Treu et al. 
2009) and the possible contribution of the environment to the 
formation of the lens (Oguri et al. 2005; Dobke et al. 2007; 
Faure et al. 2009b). To fully understand the properties of a lens 
sample, we need measure the redshift of the lenses and sources, 
probe their environments and model the mass potential of each 
lens system. 

Faure et al. (2008), presented a sample of 67 new strong 
lens candidates discovered by visual inspection of bright/early 
type galaxies with photometric redshifts < 1.0 in the COSMOS 
field. Because of the selection method, the sample po- 
tentially contains systems with exceptional individual prop- 
erties. In this paper, we study one of these lens can- 
didate systems: COSMOS 5921+0638 (RA=09 ;, 59"™21.7 S , 
DEC=+02°06'38")- The system consists of four point-like ob- 
jects that lie on top of a perfect ring around an early-type galaxy. 
According to the best-fit spectral energy distribution, a pho- 



WFPC2 F606W 



tometric redshift of z^^O^Si^, was inferred (Faure et al. 
2008). The morphology of the system suggests that it is a lensed 
AGN, the ring having been formed by the AGN's host galaxy, 
similar to that shown by RXSJ 1 131-1231 (Sluse et al. 2003) or 
PG 1 1 15+080 (Impey et al. 1998). 

We organize the paper as follows. In Sect. 2, we present 
the large imaging and spectroscopic datasets. In Sect. 3, we 
analyze the nature of the system using information from all 
these datasets. The astro me try and photometry of the objects 
in the system obtained from light profile fitting procedures 
and techniques are presented in Sect. 4. In Sect. 5, we ana- 
lyze the environment both around and along the line-of-sight to 
COSMOS 5921+0638. By using all the information gathered in 
the previous sections, in Sect. 6 we present mass models of the 
system, and in Sect. 7 we consider the different phenomena that 
can cause the observed flux ratio anomalies between the lensed 
images. Throughout the paper, we assume a flat cosmology with 
Qm=0.3 and Ho=70 km s^Mpc -1 . All magnitudes are given in 
the AB system. 



Fig. 1: HST 3" side length cut-out images of COSMOS 5921+0638. 
The left panel shows the WFPC2 F606W exposure and the right panel 
shows the ACS F814W exposure with the naming scheme selected for 
the point-like objects. 
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2. Imaging and spectroscopic dataset 

2.1. Imaging dataset from ground and space 

The COSMOS field (Scoville et al. 2007) is a square field with 
1.4 degrees side length (~2 square degrees). It was selected 
to be close to the celestial equator in a region with minimum 
extinction by dust in our galaxy (<E(b-v)>- 0.02 mag), en- 
suring maximum observability and high observation depths. It 
has been observed using multi-band imaging from space and 
Earth. Thus, we have access to data of broad wavelength cov- 
erage for COSMOS 5921+0638 and its neighborhood. Among 
these, the system was observed with the Subaru and CFHT 
ground based telescopes in the B, V, r+, i+, z+ and u* bands 
(data described in detail by Capak et al. 2007). From the Hubble 
Treasury programs (Scoville et al. 2007), we have HST/ ACS 
observations in the F814W band (the data reduction process 



Fig. 2: CFHT/Megacam and SUBARU/Suprime images of 
COSMOS 5921+0638. The 3" side length images are displayed 
from shorter (u*) to longer (z+) wavelengths from left to right, top to 
bottom. 



is described in Koekemoer et al. 2007). Additionally, we have 
access to WFPC2 F606W exposures of the system (HST pro- 
posal id: 11289, PI: Kneib); part of the Strong Lensing Legacy 
Survey (SL2S, Cabanac et al. 2007), where the system was also 
serendipitously found. In Fig. 1, we display the HST exposures 
of the system with the chosen naming scheme for the point-like 
images (A to D, clockwise from the east-most image), which 
is used hereafter. In Fig. 2 we show the CFHT and Subaru ex- 
posures. Even though the system was observed in both the ra- 
dio (1.4GHz VLA-COSMOS, Schinnerer et al. 2007) and X- 
ray (XMM/Newton-COSMOS, Brusa et al. 2007), no signal was 
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Fig. 3: Subaru/Suprime B band image of the 8'x8'field around COSMOS 5921+0638. The central lensing galaxy of the system is shown in the 
center. The galaxies with redshifts from the zCOSMOS catalog are displayed in symbol-coded redshift bins. Galaxies with measured redshifts 
from the FORS1 observations are labeled and shown with large circles. 



measured (as already mentioned in Faure et al. 2008). We sum- 
marize the properties of the imaging dataset in Table 1 . 



Table 1: Summary of the different imaging data used for this anal- 
ysis. The magnitude limits for the COSMOS observations are those 
shown in Capak et al. (2007) for the ground-based observations and in 
Koekemoer et al. (2007) for the HST observations, (a) The u* frame 
was created by combining observations taken between the years 2004 
and 2005, thus, the total exposure time of the COSMOS field in this 
band varies, (b) 27.1 is the magnitude limit for a point-like object in 
the ACS exposure and increases in brightness to 26.1 for a l"-wide 
extended object. 



Camera 


Band 


Date 


Exp. [s] 


Limit [mag] 




B 


2004-02-19 


4320 


27.3 




V 


2004-02-18 


3240 


26.6 


Suprime 


r+ 


2004-01-19 


2160 


26.8 




i+ 


2004-01-22 


2880 


26.2 




z+ 


2004-01-21 


4320 


25.2 


Megacam 


u* 


2004-2005 (fl) 


-40000 


26.4 


ACS 


F814w 


2004-04-08 


2028 


27.1 (i " 


WFPC2 


F606W 


2008-01-07 


1200 


26.5 



2.2. Spectroscopic dataset 

Galaxies in the COSMOS field have been spectroscopically 
observed as part of the zCOSMOS survey (Lilly et al. 2007). 
zCOSMOS is a large redshift survey within the COSMOS field 
using the VIMOS spectrograph installed at the VLT. The sec- 
ond data release of the survey (zCOSMOS-bright DR2, re- 
leased in October 2008; Lilly & zCOSMOS Team 2008), con- 
tains ~10,000 galaxy spectra with associated redshifts. Galaxies 
in the field around COSMOS 5921+0638 with zCOSMOS red- 
shifts (from the DR2) are displayed in Fig. 3. 

Besides the zCOSMOS spectra and redshifts, the fields of 
8 of the 67 strong lens candidates presented by Faure et al. 
(2008), including COSMOS 5921+0638, were observed with 
the FORS 1 instrument at the VLT in Multi Object Spectroscopy 
(MOS) mode as part of a follow-up program (PI: Faure, Proposal 
ID: 077.A-0473(A)). The FORS1 observations around the field 
of COSMOS 5921+0638 were obtained in April 2006 and are 
presented in this paper. They comprised a 7'x7' field centered 
on COSMOS 5921+0638 with the 1501 grism (wavelength cov- 
erage: [3300 - 6500]A, resolution: 5.54A/pixel), with a total ex- 
posure time of 1800s. The standard CCD reduction and the spec- 
tra extraction was completed using pipeline recipes provided by 
ESO 1 . The flux calibration was done using the long-slit spec- 



1 http://www.eso.org/sci/data-processing/software/pipelines/fors/fors- 
pipe-recipes .html 
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Fig. 4: Placement of the slit (shaded region) containing three (bright- 
est) images and the lensing galaxy in COSMOS 5921+0638 shown on 
top of a high resolution 3" side length ACS image cutout. The FWHM 
of the VLT FORS1 observations is shown at the top left corner of the 
image. The coarse FWHM plus the spectral dispersion direction (per- 
pendicular to the orientation of the slit) ensures that the spectrum is 
highly contaminated by the much brighter galaxy. 



trum of the LTT4816 standard star and iraf 2 routines. In Fig. 4, 
we display the position of the central MOS slit, placed on top 
of the system: its goal was to measure the redshift of the central 
galaxy, as well as obtain signal from the close north-east point- 
like objects. 



3. The nature of the system 

The lens nature of an object is confirmed with the fulfillment of 
different steps, such as the corroboration of the appropriate num- 
ber, relative brightness, and configuration of the candidate lensed 
images. The first step is, naturally, the confirmation of the line- 
of-sight alignment of the foreground "lens" and the background 
"source" with the observer by means of redshift measurements. 
However, in some cases, these can be difficult to confirm (for an 
extended discussion see Marshall et al. 2009). For instance, the 
lens or the multiple images of a source can be too faint, mak- 
ing it difficult to acquire the spectroscopy required to measure 
their redshifts and corroborate that the candidate lens and source 
are foreground and background objects, respectively. This is pre- 
cisely the case for COSMOS 5921+0638: the point-like images 
are very faint in comparison to the lensing galaxy (see Sect. 4). 
Nevertheless, an ensemble of evidence allow us to infer the grav- 
itational lens nature of the system. 

3.1. The multi-band images 

The distribution of the point-like images around the central 
galaxy shows a typical strong lens configuration: a pair of close 
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2 iraf is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 



Fig. 5: The flux-calibrated VLT/FORS1 spectrum of the central lens- 
ing galaxy towards COSMOS 5921+0638. Identified lines are labeled. 
Shaded areas denote regions of strong sky emission or absorption. 



images (A and B) and two images spread toward north (C) 
and south-west (D), similar to the configuration of the lensed 
quasars PG1 115+080 (Schechter et al. 1997), MG0414+0534 
(Hewitt et al. 1992), and WFI J2033-4723 (Morgan et al. 2004). 
Moreover, the point-like shape of the images suggests that the 
background source is an AGN. This hypothesis is also supported 
by the HST images, where a smooth structure forming a perfect 
Einstein ring below the point-like images, can be interpreted as 
emission from the host galaxy of the AGN. Based on this in- 
terpretation, we expect the point-like images to exhibit typical 
AGN emission lines such as Lya (/l=1215A) or CIV (/l=1539A). 

As shown in Fig. 2, the point-like structures have very sim- 
ilar wavelength dependence. In the CFHT/Megacam u* band 
images (limiting magnitude=26.4, central wavelength=3797A), 
where there should be very low contamination by the central 
galaxy (based on the photometric redshift z~0.5, the 4000A 
break should be at ~6000A), we are also unable to detect emis- 
sion from the point-like images. The multiple images become 
visible in the Subaru/Suprime camera B band frames (limit- 
ing magnitude=27.4, central wavelength=4459A), where their 
brightness is comparable to that of the galaxy. Toward longer 
wavelengths, their light becomes rapidly contaminated by that 
of the galaxy. This behavior suggests that the point-like images 
are related to each other. Furthermore, it provides by itself an es- 
timate of the redshift of the source, by a criterion similar to that 
present in the Lyman break technique (e.g., Giavalisco 1998) 
used to detect high-redshift star-forming galaxies. Assuming that 
our background source is an AGN, which are Lya emitters, we 
can classify it as a u* "drop-out". Thus, we can set the location 
of the Lyman break + Lya line longwards of the u* band central 
wavelength, which provides a lower limit to the redshift of the 
background AGN of z«;3. 

3.2. The FORS1 spectra 

From the slit covering a large portion of the system (see Fig. 4), 
we obtain a single spectrum. Even though it contains emission 
from both the point-like objects and the central galaxy, the ob- 
served spectrum is consistent with a typical early type galaxy 
SED at a redshift of z=0.55 1+0.001 (see Fig. 5). This value was 
derived using the Call H, Call K, G, Hy, and Mgb absorption 
lines. 

Although the slit of the FORS1 observations on top of the 
system was aligned so as to contain three of the point-like ob- 
jects (see Fig. 4), at first sight no emission line from any of them 
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Fig. 6: The COSMOS 5921+0638 central galaxy spectrum. Top: el- 
liptical template spectrum. Middle: measured spectrum of the lensing 
galaxy + point-like images. Bottom: difference between the two spec- 
tra. The thick gray lines show the continua used to scale the respective 
spectra prior to subtraction. The difference spectrum shows an evident 
residual feature at around 5050 A that we interpret as Lya- at z=3.14. 
Shaded areas denote regions of strong sky emission or absorption. 



Table 2: Galaxies with measured redshifts from the FORS1 observa- 
tions. The measured errors for the redshifts are A z=10~ 3 . 



RA 



DEC 



09:59:09.69 02:07:53.6 0.353 

09:59:12.21 02:08:25.4 0.354 

09:59:12.39 02:06:56.4 0.352 

09:59:12.41 02:07:51.2 0.352 

09:59:14.44 02:07:26.0 0.355 

09:59:15.59 02:06:12.0 0.602 

09:59:21.79 02:06:38.7 0.551 

09:59:30.27 02:04:22.0 0.687 

09:59:32.69 02:06:08.8 0.527 

09:59:33.24 02:05:47.6 0.312 



is seen. The single spectrum that could be extracted is highly 
contaminated by the much brighter galaxy in the spectral range. 
However, when subtracting a scaled elliptical-galaxy template 
spectrum (Kinney et al. 1996) from the observed spectrum (see 
Fig. 6), a small feature is observed at 5050A. We note, however, 
that the detection significance of this feature depends on the tem- 
plate selected for the subtraction. 

The asymmetric shape of this feature suggests that it is Lya 
emission. Furthermore, since the galaxy is brighter than the 
AGN, especially at longer wavelengths, and because of the pho- 
tometric u* drop-out criterion, we infer that Lya is the only 
emission line that could be seen. The Lya line is in general 
brighter than other broad emission lines and is located (in this 
case) at observed wavelengths corresponding to the fainter part 
of the galaxy's spectrum, thus, is less contaminated by the galaxy 
light. This line yields a redshift for the point-like images of 
z=3. 14+0.05. This redshift is assigned by the single line and the 
error determined by its width. 

From the spectra obtained from the remaining slits of the 
FORS1 observations, we derived the redshifts of 9 additional 
galaxies in the ~T (or ~2.7 Mpc at z=0.55) field around 
COSMOS 5921+0638. In Fig. 3 and Table 2, we report the lo- 
cation and redshift of these galaxies, 5 of which are at a similar 



redshift and relatively close to each other. This probably implies 
that they are members of a group or cluster of galaxies at za:0.35. 

3.3. Our interpretaion of the point-like objects 

If COSMOS 5921+0638 is a genuine lensed AGN, as implied 
by the typical physical configuration of the system (four im- 
ages of a source, a close pair, an Einstein ring formed by the 
host galaxy), the u* drop-out criterion suggests that the redshift 
should be z>3. The identification of a possible Lya emission 
feature using the FORS 1 spectrum is evidence that this is indeed 
the case and that the source is an AGN at a candidate redshift of 
z. s =3. 14+0.05. 

4. Astrometry and photometry of the system 

A lensing mass model depends mostly on the relative positions 
of the source images and lensing galaxy. However, other parame- 
ters obtained from a quantitative measurement of the light distri- 
bution of the objects in the system (e.g., flux ratios, ellipticities, 
position angles, see for example Faure et al. 2009a) can also be 
useful as constraints or as a comparison. It is hence important to 
accurately measure the astrometry of the objects in the systems 
and take advantage of the additional information retrieved in the 
measurement process in order to develop and interpret a reliable 
mass model of the lens. 

For this purpose, we applied the galfit software (Peng et al. 
2002) to the WFPC2 (F606W) and ACS (F814W) space-based 
observations. The software allows us to fit analytical two dimen- 
sional light profiles to the objects seen in a given CCD frame, 
in our case, the four point-like images of the background source 
and the lensing galaxy. 

4.1. The PSFs 

As an input, galfit requires the point spread function (PSF) 
of the observations to convolve it with the different analyti- 
cal profiles. We obtained the PSF model for the WFPC2 ob- 
servations using tinytim 3 . Because of our conclusion about the 
nature of the point-like images the PSF was created by as- 
suming a power law continuum with a v =-0.46 (fit from SDSS 
quasars; Vanden Berk et al. 2001) to avoid PSF-color artifacts. 
Since the WFPC2 frame of the system was created by combin- 
ing 3 dithered exposures, the modeled PSF does not take into 
account convolution effects. Nevertheless, as shown in the next 
subsection, the WFPC2 data is not used to obtain any acute lens- 
ing constraints (e.g., astrometric parameters), but only the pho- 
tometry of the objects, which is not significantly affected by con- 
volution effects. 

The ACS PSF exhibits far greater temporal and spatial varia- 
tions than the WFPC2 PSF. To obtain an accurate PSF model, 
we referred to Rhodes et al. (2006, 2007), which contains a 
statistical study that yields a measure of the focus value / 
(to a micro-metric accuracy) of the different COSMOS ACS 
exposures. Using this information for the ACS exposure of 
COSMOS 5921+0638 (/~-4.5jum) we applied the idl/tinytim 
procedures presented in Rhodes et al. (2006) 4 to generate PSF 
models, in particular choosing those in the location of the CCD 
where the system is located. The PSF was created by assuming a 
single monochromatic slice at /l=8140A (central wavelength of 
the F814Wband). 



http://www.stsci.edu/software/tinytim 
4 http://www.astro.caltech.edu/~rjm/acs/PSF/ 
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Fig. 7: Top panel: galfit de Vaucouleurs + 4 point-source fit to the ACS 
exposure of COSMOS 5921+0638: from left to right: Original image, 
fitted model, and subtraction residuals. Bottom panel: Fit of the WFPC2 
exposures with the same display scheme. 



Table 3: Parameters of the de Vaucouleurs fit to the lensing galaxy 
towards COSMOS 5921+0638. 



Filter 


Reff ["] 


e 


PA 


mag 


F606W 


0.45 (fix) 


0.148(fix) 


27.3 (fix) 


21.72+0.02 


F814W 


0.45+0.01 


0.148+0.009 


27.3+2.0 


20.310+0.009 



Table 4: Photometry of the lensed AGN images in 
COSMOS 5921+0638, magnitude values. Images labeled as in 
Fig. 1. 



Filter 


A 


B 


C 


D 


F606W 


25.60+0.08 


24.57+0.04 


24.83+0.06 


26.50+0.18 


F814W 


25.42+0.06 


24.59+0.03 


24.53+0.03 


26.19+0.10 



Table 5: Flux ratios of the AGN images calculated with the magnitude 
values shown in Table 4. 



Date 


Filter 


A/B 


C/B 


D/B 


2004/04/08 


F814W 


0.47+0.03 


1.06+0.04 


0.23+0.02 


2008/01/07 


F606W 


0.39+0.03 


0.79+0.05 


0.17+0.03 



Table 6: Astrometry obtained with galfit for the background source 
images. All quoted values are in arcseconds with respect to the lensing 
galaxy at RA=09' ! 59"""21.768 S , DEC=+02°06'38.33" 



Filter 


A 


B 


C 


D 




0.702 


0.617 


-0.067 


-0.468 


F814W . 

&DEC 


(±0.005) 


(±0.003) 


(±0.003) 


(±0.008) 


-0.132 


0.382 


0.727 


-0.507 




(±0.005) 


(±0.002) 


(±0.002) 


(±0.008) 



4.2. The light-profile fitting 

The lensing galaxy was fitted by a de Vaucouleurs profile 
(de Vaucouleurs 1948). The parameters that define this profile 
are: the effective radius (or half-light radius, R e ff), the central 
position, the ellipticity (defined as <?=^; a: semi-major axis and 
b: semi-minor axis), the position angle (PA, defined in degrees 
measured from north to east), and the magnitude. 

The four point-like images were parameterized as point 
sources, which have central positions and magnitudes as free 
parameters. Because of the lower resolution of the WFPC2 ob- 
servations, the parameters that define the de Vaucouleurs profile 
were fixed to those obtained with the ACS fitting (with the ex- 
ception of the magnitudes and central positions). In both cases, 
the background was fixed to the median value of empty regions 
selected in the field. The fitted light profiles and residuals are 
displayed in Fig. 7. The parameters obtained for the galaxy are 
summarized in Table 3 and the photometry of the point-like im- 
ages is presented in Table 4. Additionally, Table 5 shows the flux 
ratios, with respect to the bright image B, derived from these 
magnitudes. The fitted positions of the different objects in the 
system are displayed in Table 6. Since the most accurate fit was 
obtained for the ACS exposure, we quote and use throughout the 
rest of this paper the astrometry obtained with this observation. 

The error bars provided by galfit are based on the assump- 
tion that there has been a perfect fit to the data (i.e.,^ 2 =l), which 
is not the general case. Hence, to obtain error bars for the differ- 
ent parameters of the fit, we created 500 Monte Carlo realiza- 
tions for each of the HST frames and inferred the uncertainty 
from the scatter in the galfit results. These errors are in close 
agreement with those shown by a single galfit fit on the ACS 
data, but in the WFPC2 data, where the fit of the lensing galaxy 
is not as good, the errors delivered by galfit are slightly under- 
estimated. 

We attempted a light profile fitting of the ground-based ob- 
servations. However, due mainly to the large PSF in these obser- 



vations (FWHM~l") and the relatively small size of the system 
(~0.8"), the results were inconclusive. 

5. The neighborhood of COSMOS 5921+0638 

Before developing a mass model of the lens, we considered the 
influence of secondary structures along the line-of-sight to the 
source that could perturb the properties of the lensing system. 
For this purpose, we use the zCOSMOS optical group catalog 
(Knobel et al. 2009). The catalog, created using spectroscopic 
information from the zCOSMOS redshift survey, contains 800 
galaxy groups. They can be characterized by different quantities 
among which the "fudge" virial mass and the "fudge" virial ve- 
locity are the most useful for our purposes. The "fudge" quanti- 
ties are obtained using the group richness and redshift as a means 
of obtaining the virial velocity and mass, adopting the observed 
relation between these and the mock catalogs created from nu- 
merical simulations for the COSMOS field. For a detailed dis- 
cussion of the fudge virial mass we refer to Knobel et al. (2009) 
(the procedure to obtain the fudge virial velocity is analogous). 

To study the gravitational influence of the optical groups on 
COSMOS 5921+0638, we selected the 21 groups located in a 
circle of ~5' radius (or 2 Mpc at redshift z=0.551) centered on 
COSMOS 5921+0638 (see Fig. 8). The location and redshift of 
the 5 galaxies at redshift z^0.35 obtained with the FORS1 ob- 
servations (see Sect. 3.2), coincides with the most massive group 
selected (ID: 20 in the original catalog). This massive galaxy 
group at redshift z=0.354 is the only one in our selection that 
has also been identified in the X-ray regime (Finoguenov et al. 
2007). 

Using the central position and mass (or velocity dispersion) 
of the 21 groups selected, we evaluate k and y at the position of 
the lens using three different assumptions for the mass profile of 
the groups: (i) point mass, (ii) singular isothermal sphere (SIS) 
and (iii) truncated isothermal sphere (TIS). To project the k and 
y values of the different foreground and background groups (ob- 
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Table 7: Total contribution of the groups around COSMOS 5921+0638, 
assuming the different profiles. 




Fig. 8: The 21 groups from the zCOSMOS optical group catalog around 
COSMOS 5921+0638, displayed as circles. The relative sizes of the 
circles scale with the "fudge" velocity dispersion. The lensing galaxy 
is displayed in the center as a square. In the top panel, the groups are 
projected onto the observer's plane as empty circles for reference. The 
middle and bottom panels additionally show the redshift distribution in 
right ascension and declination, respectively. 



tained with each mass profile assumption) to the redshift of the 
lens, we used the methodology described in Keeton (2003) (see 
also, Momcheva et al. 2006). The point mass model assumes 
a singular total mass (which we selected as the virial mass), 
thus, there is no surface mass density outside the singularity and 
we can only observe the tidal shear y produced by this singu- 
larity. Assuming a singular isothermal sphere (SIS) profile for 
the groups, we can observe both the convergence and the shear 
(which are equivalent in this profile). Nevertheless, since the to- 
tal mass of the SIS profile diverges, the values of k derived with 
it are certainly overestimated far from the center of each group. 
Truncating the profile eliminates this issue (TIS profile). The 
virial radius of each group was chosen for this truncation. 

By adding the individual k and y values (or their vector sum 
in the case of the shear) of the 21 galaxy groups, we can de- 
termine the net influence of the environment (see Table 7). As 
expected, the total environmental convergence k produced by the 
groups, assuming a point mass potential, is zero. Assuming a SIS 
profile, we observe a total convergence /e=0.062, and for the TIS 
profile, a much smaller /e=0.010. From the latter, we can see that 





P. Mass 


SIS 


TIS 


Ktot 


0.000 


0.062 


0.010 


Jtat 


0.006 


0.005 


0.003 




14° 


-43° 


-43° 



the total k contribution is negligible. However, it is important to 
remark that because the total convergence is a sum of positive 
scalar values, incompleteness in the catalog (i.e., unidentified or 
partially identified galaxy groups) leads to an underestimation of 
this value. On the other hand, the total external shear is a vector 
sum, thus, incompleteness in the catalog may lead to an under- 
or overestimation of the value. As a result, a smaller effect on the 
total sum of the shear would be observed (even though its direc- 
tion can be compromised). Additionally, since shear is produced 
by tidal "pulls" from the mass distribution, the choice of poten- 
tial profile does not have a large influence on the final result. As 
the total shear shown in Table 7 is comparable to the same or- 
der of magnitude as the variance expected by cosmic shear, the y 
contribution exerted by the zCOSMOS groups in the catalog also 
cannot produce significant perturbations to the lens potential in 
COSMOS 5921+0638. 

6. The mass modeling of the lens 

To confirm the lens nature of the system and determine its mass 
distribution, we perform mass potential modeling of the lens by 
developing two types of mass models, grid-based and analyti- 
cal. The grid-based models have the advantage of making no 
assumptions about the mass density profile of the lens. On the 
other hand, the analytical mass model, assumes a particular type 
of mass density profile, but provides information about the lens- 
ing galaxy that can be compared to observational quantities. In 
all models, we assume that the lens plane is at z/=0.55 1 and that 
the source plane is at z s =3.14. 

6. 1 . Grid-based mass model 

The pixelens code (Saha & Williams 2004) 5 provides a grid- 
based mass model of the lens using the observed image po- 
sitions as constraints. A requirement of the minimization is 
that the arrival time order of images is known. If incor- 
rect image orders are entered, the model will either not con- 
verge or the fitted time-delay surface will determine criti- 
cal points in locations not occupied by images of the back- 
ground source. With this constraint, by trial and error, we can 
obtain the correct time delay order of the images (and the 
type of singularity on which they are placed), which in this 
case was: C(minimum)i-^A(minimum)i->B(saddle)i->D(saddle). 
Using this order, we obtain a pixelated profile for the lens. The 
time-delay contours for the correct image order are shown in Fig. 
9 (top left panel). 

We tested two models, a first model without any external 
shear and a second model with external shear. In grid-based 
modeling, the mass profile is free and a non-symmetrical mass 
distribution is inferred from the fit. However, by dividing the 
mass distributions into symmetrical and non-symmetrical com- 
ponents (as shown in Fig. 9, lower panel), we can see that the 
non-symmetric residual is rather small (for both models). 



http://www.qgd.uzh.ch/projects/pixelens/ 
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Fig. 9: Grid-based model (including external shear) of 
COSMOS 5921+0638. The top left panel shows the time-delay 
surface contours which map the location of the critical points (saddle, 
maxima and minima), and thus, the location of the images. The top 
right hand panel shows the fitted mass distribution while the bottom left 
and right panels show the symmetric and non-symmetric contributions 
to the mass distribution, respectively. Mass contours are defined as 
-2.5/ogio(M). The non-symmetric part accounts for approximately 7% 
of the total mass. 



For the model without an external shear, we obtain a mass 
profile with a slightly larger ellipticity than that of the observed 
light distribution of the galaxy. Additionally, the position angle 
of this mass distribution appears to be larger than 45° (compared 
to the 27.3° observed position angle of the light distribution). 
The amount of non-symmetrical mass in the mass profile ac- 
counts for approximately 10% of the total mass. 

The model with an external shear (shown in Fig. 9) shows an 
ellipticity similar to that of the observed light distribution with 
a position angle that is also similar to the observed one. In this 
case, the non-symmetrical residual of the mass distribution ac- 
counts for only 7% of the total mass. It is therefore plausible 
that, even though small, the contribution of an external shear is 
required to model the system. We note that, although the con- 
tribution from the asymmetrical part of the mass profile located 
south west of the center of the lens galaxy is low, this may be 
indicative of some kind of substructure within the lens. 

In both grid-based fits (with and without external shear), be- 
cause of the steepness degeneracy (e.g., Falco et al. 1985), hav- 
ing as constraints only the position of the images (lying on a 
ring) and a low ellipticity, the slope of the profile is very degener- 
ate (see Fig. 10). Nevertheless, the total mass inside the Einstein 
radius (0.71" or 4.5 kpc at the redshift of the lensing galaxy) is 
well defined and equals M# £ = 1.2xlO n Mo in the cases with and 
without external shear. 

Time-delays for both models are relatively similar but have 
large error bars. They are of the order of half a week for the 
longest delay (C to D) and of the order of a couple of hours for 
the shortest (A to B). 

Grid-based modeling (Saha & Williams 2004) provides in- 
sights on the mass distribution of the lens and a familiariza- 
tion with the system. However, being grid-based, it naturally has 
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Fig. 10: Radially enclosed mass obtained with the grid-based model that 
included external shear. The segmented lines display the effective radius 
fitted for the lensing galaxy and the measured Einstein radius of the 
system, respectively. 



Table 8: SIE fit to the lensing system. All observational constraints 
were fixed to the observed values. The galaxy's central position was 
allowed to move within 0.01", which is the estimated uncertainty. 





SIE 


SIE+y 


Ax ["] 


0.01 


-0.007+0.002 


Ay ["] 


0.005 


-0.005+0.002 


<r [km/s] 


190.7 


189.3 ± 0.4 


e 


0.15 (fixed) 


0.15 (fixed) 


PA 


27.3° (fixed) 


27.3° (fixed) 


r 




0.038 ± 0.002 






107.4 ± 1.3 


xt 


-100 


0.3 



enormous freedom in defining the mass profile and does not al- 
low a quantitative comparison with observations. The final pixe- 
lated maps reflect the highest probability model from the ensem- 
ble of models tested (100 in our case). 

6.2. Analytical mass model 

To obtain quantitative information about the mass model param- 
eters and compare this with the information obtained by pixe- 
lens, we attempt to develop an analytical model of the system 
using lenstool 6 (Kneib 1993). lenstool does source and image 
plane parameter based x 1 minimization using a Bayesian algo- 
rithm (prior based minimization) and a Markov Chain Monte- 
Carlo (MCMC) process that samples the probability distribution 
by random variation of the parameters. This parameter space 
scanning, prevents the minimization being trapped in local x 1 
minima and allows a very robust result with minor speed trade- 
off. The details of the algorithm and the underlying statistics are 
described in lullo et al. (2007). For our mass models, we apply 
the image plane minimization algorithm and use the observed 
positions of the AGN images as constraints, with a conservative 
error of A=0.01". This error accounts for imperfections in the 
lensing galaxy's position fit because of the overestimation of the 
flux in its central cusp. 

Constraining the galaxy's central position, ellipticity, and 
position angle to be the observed values and fitting a singular 



6 http://www.oamp.fr/cosmology/lenstool/ 
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Fig. 11: lenstool's SIE+y fit to the lensing galaxy. Observed (asterisks) 
and resulting (triangles) image positions from the mapping of the posi- 
tion of the fitted background source (circle). The dashed-dotted curve 
shows the lens plane critical curve and the central continous diamond 
shaped curve shows the source plane caustics. The long-dashed and 
continous ellipses denote the observed and fitted lensing galaxy posi- 
tions, respectively. The short-dashed line shows the orientation of the 
external shear required by the model. The bottom panel shows a zoom 
to the center of the system. 



isothermal ellipsoid (SIE) yields a poor fit to the data (j^^lOO), 
as expected from the results of the grid-based modeling. 

When we add an external shear to our SIE model, we obtain 
far more accurate results (as displayed in Table 8 and Fig. 11). 
Figure 1 1 shows the comparison between the observed AGN im- 
age locations (asterisks) and the "best fit" image locations (trian- 
gles). These "best fit" image locations come from the mapping 
of the average location of the minimized source position. The 
zoom in the bottom panel, allows us to see the position of the 
source with respect to the caustics produced by the mass poten- 
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0.05 0.1 0.15 0.2 0.25 
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Fig. 12: Degeneracy between the ellipticity of the lensing galaxy and the 
external shear. The contours show the confidence regions (68%, 95% 
and 99%) for the fitted parameters. 



tial. The total mass within the Einstein radius matches the value 
inferred by the grid-based mass model (Mr e ^1.2 x 10 n M o ). 

As seen in Table 8, the model needs to include a small ex- 
ternal shear (y=0.038) to reproduce the observed image con- 
figuration. The small strength of this shear is compatible with 
us not observing a strong influence from the environment. 
Nevertheless, its strength is somewhat larger than that implied 
by the group analysis. This may be indicative of: i) a misalign- 
ment between the light and the mass distribution of the lens- 
ing galaxy, ii) a contribution from individual galaxies close to 
the main lensing galaxy, and/or iii) a more complex descrip- 
tion of the lens potential (e.g, non-isothermal slope or a mis- 
aligned dark matter halo). A misalignment between the light and 
mass distribution in this case is not a likely explanation since 
the model does not provide a good fit even when the ellipticity 
or position angle of the mass distribution are allowed to vary 
(±10° with respect to the position angle of the light distribu- 
tion; see Keeton et al. 1997). We also examined the degeneracy 
between the ellipticity of the mass distribution and the exter- 
nal shear. To achieve this, we scanned the x 1 surface produced 
when these parameters are allowed to vary (in addition to the 
velocity dispersion, and central position), while the position an- 
gles of the mass profile and external shear are fixed to their ob- 
served and best-fit model values, respectively. Figure 12 illus- 
trates the clear degeneracy between these two parameters, as the 
necessary strength of the external shear decreases as the elliptic- 
ity of the mass profile decreases. However, the best-fit model 
value for the ellipticity (es;0.13) remains consistent with that 
measured from the light profile of the lensing galaxy (e=0. 148), 
and a shear of strength y>0.01 is required at 99% confidence. 
Regarding the contribution from individual galaxies, at the time 
being, there is no spectroscopic measurement of the immediate 
neighbors. Nevertheless, the three brightest and closest neigh- 
boring galaxies (<6" in projection to the observational plane) 
are aligned in an ~ 105° east from north direction, in close agree- 
ment with the modeled external shear. Their photometric red- 
shifts indicate that they are not within the same plane as the 
main lensing galaxy (zi=l.l, Z2=0.8 and Z3=1.2). Assuming sin- 
gular isothermal sphere profiles for these three galaxies and pro- 
jecting their contribution (from their respective photometric red- 
shifts) onto the lens plane (at z/=0.551), they could reproduce 
the required external shear if they each have velocity dispersions 
of <x~300 km s . Even when velocity dispersions of that or- 
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Table 9: Convergence (k), shear (y), and magnification (jj) shown by 
the potential (SIE+y) at the location of each lens plane image. The last 
column shows the flux ratio of each image with respect to image B 
inferred for the magnification values fj. 



Image 


K 


r 




Fl. Ratio (JB) 


A 


0.472 


0.508 


47 


0.70 


B 


0.509 


0.506 


67 


1.00 


C 


0.498 


0.469 


31 


0.46 


D 


0.548 


0.525 


14 


0.21 



der have been measured in galaxies, the low brightness of these 
candidates (m/ ; i=24.1, m/2=23.2 and W/ j 3=24.5), makes it an 
unlikely possibility in these cases. 

By considering an analytical model, different properties of 
the system can be checked both in the lens and the source plane. 
The fit of the potential can be used to measure the values of con- 
vergence (k), shear (y), and, thus, the magnification (//) expected 
at different positions in the lens plane. With these modeled mag- 
nification values at the position of the images, we can then infer 
the expected flux ratios for the different images. The k, y, and 
values and the flux ratios (with respect to image B) of the differ- 
ent images are displayed in Table 9. Comparing these modeled 
flux ratios with the observed ones (see Table 5), we can see a 
discrepancy. The possible explanations of this discrepancy are 
discussed in Sect. 7. 

The magnification values displayed in Table 9 can also 
be coupled with the measured photometry of the individual 
point-like images (Table 4) to obtain an estimate of the intrin- 
sic brightness of the background source. With this, we obtain 
an approximate absolute magnitude of the source of: — 
17.5 and M%uw — 18-0. Even though these values are the me- 
dian of the absolute magnitudes (with no k-correction) obtained 
for each image, which are affected by uncertainties related to 
flux ratio anomalies and the assumption of a SIE mass pro- 
file, they reflect the low luminosity nature of the source AGN 
(LLAGN; e.g., Storchi-Bergmann et al. 1995; Bower et al. 1996; 
Cid Fernandes et al. 2004). 

6.3. Mass-to-light ratio 

As shown previously, theoretical mass modeling (both grid 
based and analytical) can provide an accurate estimate of the 
total mass enclosed within the Einstein radius. Using aperture 
photometry within this Einstein radius (i.e., #£=0.71"), we also 
measured the amount of light enclosed within it, to compare it 
with the enclosed mass. This measurement was performed on 
the fitted galfit de Vaucouleurs profile, to avoid the contribution 
from lensed images. 

We measured an ACS F814W aperture magnitude of 20.88 
for the lensing galaxy. Using a k-correction of 0.45 for an ellip- 
tical galaxy at redshift z=0.55 in the F814W filter and an abso- 
lute solar magnitude of M 0Fgl4W =4.53 7 , we obtain Lfsuw~4.3 x 
10 10 L Gf . 8]4lv . This, together with the enclosed mass (1.2 x 
1O U M , see Sect. 6.1 and 6.2), leads to a mass-to-light ratio 
of M/L F guw~2-7Me/Le muw . The WFPC2 F606W observations 
yield an aperture magnitude of 22.28 within the Einstein radius. 
Using a k-correction of 1.33 and an absolute solar magnitude of 
M o F6 mw=4-74, we obtain L F60 6w~3.3 x 1O 1O L 0RSO61V , leading to a 
mass-to-light ratio of M/L F60 6w~ 3.6M G /L 0rao6lv . 



7 obtained from http://www.ucolick.org/~cnaw/sun.html 



Fukugita et al. (1995) show typical rest-frame elliptical 
galaxy colors of B-F814W*2.23 and B-F606W-1.29. We 
couple these colors with our aperture photometry (F814W 
and F606W) data and an absolute solar magnitude of 
M#=5.33 to retrieve an estimate of the B-band luminosities 
for the lensing galaxy of Lb fsuw ~1.19 x 1O 1O L 0B using the 
F814W data, and L BfW6H ,«1.75 x 1O 1O L 0B using the F606W 
data). These lead, respectively, to: M/Lb fsuw ~1OAM q /L 0b 
and M I Lb F606W ~6.9MqI L &B which yield an average value of 
M/Lb~8.5+1.6. These values agree with those measured by 
Keeton et al. (1998) for the mass-to-light ratio within the 
Einstein radius of 10 elliptical lens galaxies (M/Lg~40 for 
galaxies at redshift z~0.55). 



7. Anomalous flux ratios 

From the observed flux ratios shown in Table 5, we can see that 
the ratios of the brightness of image A to B, and of D to B do 
not vary significantly (within the uncertainties) among the differ- 
ent filters and the different epochs (A/B 2004 JF8 14=0.47+0.03 
and A/B 2008_F606=0.39+0.03; D/B 2004_F814=0.23+0.02 
and D/B 2008_F606=0. 17+0.03). However, the C/B flux ratios 
show variations between different filters and/or over time (C/B 
2004.F8 14=1.06+0.04 and C/B 2008_F606=0.79+0.05). 

Additionally, comparing the theoretical flux ratio values 
(Table 9: A/B=0.70, C/B=0.46, and D/B=0.21) with the ob- 
served flux ratios (Table 5), a disagreement can be seen. Even 
more generally, any symmetric model that we choose for the lens 
potential, expects images A and B to be the brightest images. 
They are the merging images of the lensing system and, conse- 
quently, the images with the highest magnifications. However, 
this is not the observed case, since image C is as bright as image 
B and brighter than image A in the ACS observations. 

Several phenomena can explain the flux differences be- 
tween the lensed images from different bands and different 
epochs and their comparison with theoretical expectations, such 
as: galactic extinction, substructure, intrinsic variability of the 
background source, or microlensing (e.g., Anguita et al. 2008a; 
Yonehara et al. 2008): 

- Galactic extinction: Galactic extinction is a static phe- 
nomenon (on human timescales) in which dust in the lensing 
galaxy absorbs the blue light originating in the background 
source images. This implies that lensed images located be- 
hind different column densities of dust appear to have color 
perturbations, in particular, those behind a larger dust col- 
umn densities would appear redder. For this particular sys- 
tem, we obtain F606W-F814W colors of: 0.18, -0.02, 0.30, 
and 0.31 for images A, B, C, and D, respectively. In terms 
of the galactic extinction interpretation, there would be less 
dust accumulated over the projection of the lensing galaxy in 
the location of image B than in the others. Assuming galac- 
tic extinction (Cardelli et al. 1989) with R v =3.1, a differen- 
tial column density of A/V//~l .8 x 10 21 would be required to 
explain the color difference between the anomalous images 
C and B (Af8i4w-,F606w=0.32). 

- Mass substructures: In a similar way to galactic extinc- 
tion, mass substructures can produce static anomalies in 
the measured flux ratios. A galaxy with mass substruc- 
tures (of the order of 10 6 ~10 10 M Q ), i.e., non-symmetrical 
luminous or dark matter components, or even unseen satel- 
lite galaxies, might modify the lensing potential, particu- 
larly the flux ratios (Kochanek & Dalai 2004). In the case 
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of COSMOS 5921+0638, even though it is small, the grid- 
based model of the lens shows a non-symmetric part of the 
distribution, which might be a hint of substructure in the lens. 
This phenomenon could explain, for example, the anoma- 
lously bright image C. However, because of the chromatic 
variation seen in the flux ratios, it cannot explain by itself the 
observed flux anomalies. If substructures and galactic extinc- 
tion are responsible for the flux ratio anomalies, additional 
observations of the system should not display brightness and 
color variations with respect to the present data set. 

- Intrinsic variability: AGN vary their brightness over periods 
of years, months, weeks, and even days, usually becoming 
bluer as they become brighter. Unfortunately, the HST ob- 
servations, although taken in different filters, are also taken 
at different epochs, which causes degeneracies in the inter- 
pretation of the chromatic effect. As the AGN images are 
separated between each other by time-delays (of the orders 
of hours to half a week, according to the mass modeling), 
the intrinsic variation of the AGN should be on timescales 
shorter than these to be able to see the anomalously bright 
image C, brighter than the close pair A and B. This is un- 
likely because quasars (statistically) require timescales of the 
order of years to produce brightness variations of ~0. 1 mag- 
nitudes (for more details, see the structure function calcu- 
lations of Vanden Berk et al. 2004). Even if short timescale 
variability was occurring in this system, these rapid bright- 
ness fluctuations imply that it is unlikely after a period of 
4 years, that slightly different flux ratios, although the same 
overall brightness order, are observed. 

- Microlensing: Microlensing is also a temporal phenomenon 
in "macro"-lensed AGNs. In this case, stars in the lens- 
ing galaxy produce an additional (de)magnification of the 
"macro"-lensed images of the source. Because of the pro- 
jected location of the background source images in the lens- 
ing galaxy, microlensing is expected in this system. Much 
like intrinsic variability, (de)magnification by microlensing, 
in the general case, produces a blue excess: the blue emis- 
sion from the background source originates in a smaller re- 
gion, therefore, it is more effectively affected than the outer 
regions. As before, since the observations in different bands 
are taken at different epochs, no strong conclusions can be 
made from the chromatic point of view. The main differ- 
ence with intrinsic variability /time-delay induced anomalies, 
is that in microlensing, there is no correlation between the 
multiple AGN images. Additionally, the projected velocities 
of the stars in the lensing galaxy, the lensing galaxy itself, 
and the background source set the microlensing timescales 
to be of the order of months to decades, compatible with the 
observed variations. For the microlensing interpretation, and 
considering that the measured flux in image C is higher than 
expected (with respect to image B), the flux anomaly that 
we observe is caused either by: (i) microlensing magnifica- 
tion of the flux from image C, mainly during the observa- 
tions from 2004, (ii) de-magnification of image B in 2004, 
(iii) magnification of image B during the 2008 observations, 
or (iv) a combination of all phenomena. As image B is lo- 
cated on a saddle point (i.e., k + y > 1), its de-magnification 
due to microlensing is more likely than for the other images 
(Schechter & Wambsganss 2002), which supports the sec- 
ond interpretation. In any case, since the four lensed images 
are located in regions of similar optical depth, variations in 
the brightness of all images should be seen in time. 



8. Summary and conclusions 

We have studied in detail the gravitationally lensed system 
candidate COSMOS 5921+0638 discovered in the COSMOS 
field (Faure et al. 2008). Using data derived from the COSMOS 
observational dataset plus additional HST-WPFC2 and VLT- 
FORS 1 observations, we have obtained spectroscopic, astromet- 
ric, photometric, and morphological parameters of the objects in 
the field. The FORS1 observations allowed us to measure the 
redshifts of 10 galaxies in the field of COSMOS 5921+0638, 
including the central galaxy of the system at z=0.551 ±0.001. 
Coupling different criteria with the observations available, we 
conclude that this galaxy is lensing a background AGN and its 
host galaxy located at a tentative redshift z=3. 14+0.05. 

We have estimated the lensing contribution of the environ- 
ment to the lensing system using the zCOSMOS optical group 
catalog. We have assigned different spherical mass profiles to 
the groups, constrained by their virial mass estimates. The anal- 
ysis yields as a result a low environmental contribution at the 
lens position (K groups ~Q.Ql and yg, oli/ „~0.005). It is important to 
note that this analysis strongly depends on the completeness of 
the optical group catalog because the contribution from single 
undetected groups may influence the results. Furthermore, the 
individual group mass estimates in the catalog may be underes- 
timated if the full mass structure of individual groups has not 
yet been identified. The forthcoming release of zCOSMOS 20K 
sample will allow us to address these issues in detail. 

We have investigated the mass profile of the lensing system 
using strong lens modeling, both grid-based and analytical. Both 
types of modeling require a small but non-negligible external 
shear (y=0. 038+0. 002) to reproduce the images' configuration. 
This presence of external shear cannot be fully explained by per- 
turbations caused by either galaxy groups in the zCOSMOS op- 
tical group catalog or by realistic M/L ratio galaxies in the sur- 
roundings of the system. The requirement of this external shear 
may therefore be indicative of incompleteness in the environ- 
mental information we have at hand (i.e., incomplete identifica- 
tion of the structure of the groups), the need to consider more 
complex mass potentials (e.g., a non-isothermal slope or a mis- 
aligned dark matter halo), a contribution from undetected indi- 
vidual galaxies, or a combination of any of these possibilities. In 
a forthcoming paper (Faure et al. 2009), we investigate the dark 
matter halo alignment with the light profile for the rest of the 
COSMOS sample. In any case, for COSMOS 5921+0638, both 
measured and modeled external contributions to the potential are 
quite small; thus, the mass of the galaxy inside the Einstein ra- 
dius has a low environmental contamination, making this system 
an interesting candidate for follow-up high resolution spectro- 
scopical observations to measure the kinematics of the galaxy 
inside the Einstein radius. 

The lens modeling of COSMOS 5921+0638 shows time- 
delays of the order of hours to days between the different im- 
ages and a total magnification of yU^150. Comparing the ob- 
served brightness of the images with the magnification in- 
duced by the fitted mass profile infers the low luminosity na- 
ture of the background AGN with an intrinsic brightness of 
M — 17.5. Additionally, by comparing the lensing mass with the 
observed light enclosed within the Einstein radius (~4.5 kpc), 
we have inferred a mass-to-light ratio of the lensing galaxy of 
M/L B «8.5+1.6. 

Flux anomalies are observed between the different lensed 
images. With the available dataset, it is impossible to determine 
a unique reason for these variations. This is mainly because the 
observations available that can be used to measure accurate pho- 
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tometry for the point-like images are in different bands and at 
different epochs. This does not allow us to ascertain whether the 
phenomenon is dynamic or static. With this in mind, our pre- 
ferred explanations of these incompatibilities are microlensing 
and/or substructures coupled with differential dust extinction in 
the lens, with microlensing in general being the most natural ex- 
planation of flux ratio anomalies in the optical range. 

By completing follow-up observations of the system, either 
the phenomenon or phenomena responsible for these anomalies 
could be found and quantified. Multi-band photometry of the 
system at a single epoch would allow us to measure the color of 
the different images, free of possible temporal anomalies. In ad- 
dition, by repeated multi-band observations at different epochs, 
the possible temporal flux variations could be proven and quan- 
tified. If active microlensing is present in this system, the addi- 
tional resolution power brought by microlensing could infer in- 
teresting properties about the background source. If the anomaly 
is static, a measure of possible substructures and projected dust 
densities in this lens can be accomplished. 

The lensing-induced magnification of the faint high-redshift 
background AGN in COSMOS 5921+0638 allows us to study 
the properties of an object that would have otherwise remained 
undetected in the COSMOS field. Even though low luminosity 
AGN (LLAGN) are expected to be far more common than high 
luminosity AGN, only a few of these objects have been observed 
beyond the local universe. 

The small separations and the low brightness of the images 
of the background AGN in COSMOS 5921+0638 make it an 
ideal candidate for additional space-based and/or new genera- 
tion optics ground-based follow-up observations to address the 
remaining open issues and understand the nature of the back- 
ground AGN. 
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